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Nonlinear optical methods are becoming ubiquitous in many areas of modern photonics. They 
are, however, often limited to a certain range of input parameters, such as pulse energy and average 
power, since restrictions arise from, for example, parasitic nonlinear effects, damage problems and 
geometrical considerations. Here, we show that many nonlinear optics phenomena in gaseous media 
are scale-invariant if spatial coordinates, gas density and laser pulse energy are scaled appropriately. 
We develop a general scaling model for (3+l)-dimensional wave equations, demonstrating the invari¬ 
ant scaling of nonlinear pulse propagation in gases. Our model is numerically applied to high-order 
harmonic generation and filamentation as well as experimentally verified using the example of pulse 
post-compression via filamentation. Our results provide a simple recipe for up-or downscaling of 
nonlinear processes in gases with numerous applications in many areas of science. 


Nonlinear interactions of intense short laser pulses with gaseous media form the basis behind a wealth of inter¬ 
esting phenomena such as multiphoton ionization and plasma formation [^, spectral broadening (which can be 
used for pulse compression [3]-[^), harmonic generation and wave-mixing [^, as well as the creation of attosecond 
pulses and the formation of electron or ion beams [^. Advances in femtosecond laser technology constantly 
yield shorter pulses, higher pulse energies, and higher repetition rates [9}|TT]. However, to fully explore this newly 
available parameter regime, which gives access to e.g. faster time scales and higher intensities, is often challenging 
because of damage problems, additional (unwanted) nonlinear effects, or geometrical restrictions. We illustrate 
this challenge for two important applications of nonlinear optics, filamentation in gases used e.g. for laser pulse 
compression, and high-order harmonic generation (HHG) providing the basis for attosecond science. 

The propagation of an intense short laser pulse in a transparent medium induces nonlinear effects caused e.g. 
by the intensity dependence of the refractive index. When self-focusing due to the Kerr effect balances defocussing 
caused by diffraction and plasma generation, a filament can be created. In addition, self-phase modulation and self 
compression may take place in the filament, resulting, possibly after further compression, in ultrashort pulses close to 
the fundamental limit of a single cycle [^. Forming a filament requires a certain power, known as the critical power 
for self-focusing (Uni. At slightly higher power, limitations arise and multiple filaments are created [^. Different 
attempts were suggested to increase the output energy ^2 16-^. However, pulse compression using filaments (or 
similarly hollow fibers) is still limited to pulse energies of typically a few mJ [^[^ , which is approximately two to 
three orders of magnitude below the maximum pulse energies available from today’s femtosecond laser sources. To 
scale up pulse post-compression into the 100 mJ-regime and above will enable the production of few-cycle pulses 
with unprecedented peak power, opening the door to new applications. 

Our second example pertains to high-order harmonic generation, which occurs when intense short laser pulses 
interact with a gas of atoms or molecules at an intensity of ^ 10^^ W/cm^ [^. This process leads to the formation 
of attosecond light pulses, which can be used for pump-probe studies of ultrafast electron dynamics [^. A major 
limitation of attosecond science is the low photon flux available [^. Since the early days, a strong effort has been 
devoted to optimize and scale up HHG [25}|^ aiming for an efficient conversion of high laser pulse energies into the 
extreme ultraviolet (XUV). In spite of this effort, propagation effects and geometrical considerations have limited 
the useful input laser pulse energy and only a few groups have employed pulse energies exceeding 10 mJ [28}|^ . In 
the opposite direction, progress in laser technology now enables the generation of laser pulses with /rJ energies at 
MHz repetition rates [^. In this regime, macroscopic phase matching issues have limited the conversion efficiency 
into the XUV and only recent attempts point towards a solution of this problem [M|[^ . An efficient down-scaling 
of HHG will enable compact attosecond sources at high repetition rate, of great interest for studies of time-resolved 
photoelectron emission processes on surfaces and for coincidence spectroscopy in atoms and molecules. 

Taken together, these examples illustrate the strong need for a general methodology that enables up- or down- 
scaling of nonlinear processes in gases. Here we present such a methodology and introduce a set of general scaling 
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Figure 1: Scaling principle: Illustration of scale-invariant nonlinear optics. A laser pulse is focused into a gas 
medium with length L and density p. Nonlinear propagation effects lead to a modification of the spatiotemporal 
pulse profile (bottom). Identical spatiotemporal modifications can be expected if a more intense laser pulse is 
focused weaker (to reach the same intensity) into a larger medium with length rfL and lower density p/r]‘^ (top). 
Note that in the up-scaled configuration, the far-held beam diameter is smaller and the amplitude larger at identical 
distance 2 ; from the focus. 


relations directly derived from basic propagation equations for ultrashort laser pulses. We verihed our theoretical 
predictions by performing numerical propagation calculations for the two phenomena described above, hlamentation 
and HHG. We show how these processes can be invariantly scaled to laser pulse energies well above the 100 mJ 
level, with no fundamental upper limit. Moreover, we experimentally verihed the invariant scalability of pulse com¬ 
pression via hlamentation within a driving laser pulse energy range exceeding one order of magnitude. Our scaling 
formalism is simple and general, and opens up completely new parameter regimes for nonlinear optics in gaseous 
media and more generally for ultrafast science. 


Scaling principles 

We illustrate our scale-invariant nonlinear optics framework using general wave equations. Nonlinear pulse prop¬ 
agation in gases (including generation of new frequencies) is usually treated using wave equations in scalar and 
paraxial approximation, which can be directly derived from Maxwell’s equations. Such wave equations describe 
electromagnetic waves propagating in one direction, exhibiting only small angles relative to the optical axis. With¬ 
out any limitation of the spectral bandwidth and thus of the minimum pulse duration, the propagation equation 
for the electrical held in frequency representation E{r^z^uj) = exp(icL;t)F^(r, 2 ;, t)dt, usually referred to as the 
Forward Maxwell Equation [^, can be written as: 


dz 
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Here, k{uj^p) = n{uj^p)ujlc denotes the wave number with angular frequency cj, refractive index n = n{uj,p), and 
speed of light in vacuum c. p is the gas density, Fnl is the frequency representation of the nonlinear polarization 
induced by the electric held E and eg is the vacuum permittivity. For short pulse propagation, exact knowledge 
of the refractive index n, e.g. in the form of a Sellmeier equation is required. For pulse propagation in the visible 
and near infrared spectral region, k{uj,p) is real-valued, but linear absorption can easily be included by a complex 
wave number. For the sake of simplicity, we consider linear polarization and rotational symmetry and thus a single 
radial coordinate r, although our theory does not require these simplihcations. The transverse Laplace-operator in 
equation 0 then becomes jdr^ + 1/r • djdr. Via the nonlinear polarization a large number of nonlinear 

interactions can be considered, such as self-focusing, self-phase modulation, held ionization, harmonic generation, 
and plasma-defocussing. 

For propagation in vacuum, the right hand side of equation 0 vanishes and k(uj^p) k{uj^0) = uj/c. We now 
introduce the held £ = Eex.p[—iuj z /c] and rewrite equation 0 : 


y 

dz 



£ = 0 . 


( 2 ) 


The change of fields from E to £ formally corresponds to a transformation of equation 0 from the laboratory frame 
to a frame moving at the vacuum speed of light c ~ 


37 . It should be noted that £ is an electric held, not an envelope. 
No envelope approximations and thus no restrictions on the spectral bandwidth are made. Equation 0 is invariant 
under the following transformations: r —^ pr and 2 : — p‘^z (see Table 1), where 7 ^ is a scaling parameter. If f (r, z) 
is a solution to the wave equation, £{rjp^zjr^) is a solution as well. Eor monochromatic waves, one prominent 
solution of equation 0 is the Gaussian beam. The scaling is obvious for the characteristic spatial parameters of 
the Gaussian beam, i.e. the beam radius IFo and the Rayleigh length zr: Wq —^ pWo and zr —While the 
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Table 1: Scaling relations: This table lists the scaling relations derived in this work, p^v and Zcr denote the critical 
power and the distance, at which an initially collimated beam collapses due to self-focusing. Sq and Tg denote the 
harmonic pulse energy and the conversion efficiency into harmonic order q. 


Gaussian beam is just one possible solution to equation Q, more generally, any kind of beam that can be described 
by this wave equation is scale-invariant under the above specified transformation. 

These basic scaling principles can be generalized to ultra-short laser pulse propagation in gases and a wide 
range of nonlinear interactions, if the medium density and the input laser pulse energy Cin are included as scaling 
parameters. By introducing £ and Pnl = exp[—icjz/c] into equation W, we obtain: 
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where p) = k{LO, p) — k{u, 0) is proportional to p and describes pulse dispersion [see supplementary information 
(SI) for details]. By neglecting the weak pressure dependence of k{u,p) in the denominator of the diffraction term, 
the left hand side of equation @ is invariant under the above transverse and longitudinal scaling transformations, 
if simultaneously the gas density is scaled, i.e. p —Similarly, the nonlinear polarization and consequently 
the right hand side of equation ^ is proportional to gas pressure p for a wide range of nonlinear interactions 
(throughout the manuscript, we assume p oc p, taking into account a constant temperature). Finally, the input 
energy 5in, proportional to the radial (and temporal) integral of the absolute square of the input field, needs to be 
scaled as 5in —to ensure that the field amplitude, which affects Pnl is kept constant under the scaling 
transformation. The output pulse energy ^out, proportional to the integral of the absolute square of the field at the 
end of the medium, follows the same scaling: Sout —^ This scaling applies as well to the generation of new 

frequencies, as shown for the case of HHG below. In practice, the geometrical scaling can be achieved by changing 
the focusing geometry (focal length and/or beam diameter before focusing) as well as the medium length. It should 
be noted that the transformation to the moving frame, leading to equation (§, was performed to illustrate the 
scaling principles, but does not constitute a general limitation of the formalism. The scaling itself is independent 
from the reference frame. 

According to the above relations (see also Table 1), any spatiotemporal modifications of the field induced by 
diffraction, dispersion, or a nonlinear process that is proportional to pressure, are scale-invariant. In practice, an 
optical process in a gas medium, defined by a nonlinear effect and certain input parameters (pulse energy, gas 
pressure, focusing geometry), can be up- or down-scaled to different pulse energies without changing its general 
characteristics. Furthermore, our scaling formalism preserves the carrier-envelope phase (CEP), which only changes 
because of linear and nonlinear (e.g. self phase modulation) propagation effects, both of which are scale invariant. 
This implies that strongly CEP-dependent processes such as single attosecond pulse generation can be invariantly 
scaled. The scaling principle is illustrated in Eig. using the example of beam reshaping under the influence of 
nonlinear propagation and applied below to filamentation and attosecond pulse generation. 


Filamentation 

A prominent example where several nonlinear propagation effects play a critical role, is filamentation [^. A 
characteristic parameter influencing the onset of filamentation dynamics is the so-called critical power for self- 
focusing Per = Ncr^?/ 47rnon2. Here, N^r is a constant depending on the spatial beam shape {N^r = 1.896 for a 
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Figure 2: Scaling filament at ion and HHG: a, b, Simulated spatiotemporal intensity distributions (normalized 
individually) in a focused laser beam in Ar for three different positions along the propagation axis and two input 
parameter sets, scaled according to the presented scaling relations (a: r = 20 fs, ein = 2 mJ, p = 1.2 bar, Wq = 40 /rm, 
b: r = 20fs, 5in = 128mJ, p = 18.75mbar, Wq = 320/rm). e, f. Simulated spatiotemporal intensity distributions 
for high-harmonic emission (above 31.5 eV) in Ar at three positions within the nonlinear medium (e: r = 10 fs. 
Sin = 62.5/iJ, p = 256 mbar, Wq = 10.6/rm, L = 2 mm, f: r = lOfs, = 16 mJ, p = 1 mbar, Wq = 169.6/rm, 
L = 0.51 m). For both filamentation and HHG, the longitudinal position is specified with respect to the position of 
the geometrical focus; in a and b in units of the respective Rayleigh lengths and in e and f in units of the length of 
the generation medium L. c, characteristic length, i.e. filament and gas cell length, respectively (blue, left axis) and 
gas pressure (red, right axis) as a function of rj and ein- t] was arbitrarily set to unity for s-m = 1 mJ. d. Integrated 
relative scaling error for the filament scaling presented in a and b for intensity (dots) and fluence (circles) as defined 
in Methods. 


Gaussian transverse profile [^), A is the laser wavelength, no the refractive index at the central frequency and n 2 
the nonlinear refractive index. Since n 2 is to very good approximation, proportional to the gas density, Pcr ^ V^Pcr^ 
thus following the same scaling relation as 5in, i.e. the critical power increases linearly with laser pulse energy. It 
can also be shown that the distance Zcr at which an initially collimated laser beam collapses due to self-focusing, 
scales quadratically with the initial beam size (i.e. ^ confirming that the scaling transformations 

remain valid under nonlinear propagation conditions. 

We performed a more rigorous verification of our scaling model by numerically simulating filamentation with a 
state-of-the-art pulse propagation code (see Methods for details). Figure [^(a) and (b) illustrate filamentation in Ar, 
using a 20 fs input pulse centered at 800 nm and two different parameter sets, where parameter set (b) corresponds 
to the up-scaled parameters (77 = 8 ) of parameter set (a). In Fig. |^(a) and (b), the spatiotemporal intensity 
distribution is shown for three positions along the optical axis. In both cases, typical filamentation characteristics 
like conical emission and temporal self-compression can be observed. Despite the very different pulse energies 
[ 77 ^ = 64 times larger for parameter set (b)] and transverse scales (77 = 8 times larger), only minor differences are 
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visible demonstrating the validity of the scaling model for filament at ion. Figure [^(c) illustrates how experimental 
parameters like input energy, gas pressure and filament length (defined here as the propagation length over which 
the intensity on the optical axis exceeds 5 • 10^^ W/cm^) scale with 77 . 

Figure |^(d) shows a numerically extracted relative scaling error, representing the deviation from perfect seal- 
ability for output intensity (dots) and fluence (circles) as a function of e-m- For each pulse energy, the error was 
calculated by comparing the output intensity (or fluence) to that obtained with four times larger pulse energy (see 
Methods for details). While the scaling error is negligibly small for pulse energies well above 1 mJ, thus indicating 
no fundamental upper scaling limit, a clear deviation from perfect scaling appears for small pulse energies. These 
deviations can be mainly attributed to avalanche ionization (see SI). 

Attosecond Pulse Generation 

As a second example, we consider HHG in gases, which leads to the formation of attosecond pulse trains or single 
attosecond pulses. Although both cases are easily encompassed in our scaling model, for simplicity we concentrate on 
pulse trains. HHG in an extended nonlinear medium can be described in two steps: first, the laser pulse propagates 
through the nonlinear medium, inducing a polarization Pq = 2dqp^ at multiple, odd-order harmonic frequencies, 
where dq is the single atom nonlinear dipole moment. Second, the harmonic field Eq is generated from the induced 
polarization. The propagation of £q = Eqex.-p{—iujz/c), where uj now denotes the harmonic frequency, can be 
described by equation with ^nl being replaced by Vq. Since both the fundamental and the harmonic fields 
follow scale-invariant propagation equations, HHG is invariant under the scaling transformations. Gonsequently the 
harmonic output pulse energy Cq follows the same scaling Cq —> This implies that the conversion efficiency 

Tg = Cqjsin is scale-invariant. In other words, the same conversion efficiency can be expected for HHG driven by 
intense laser pulses, with a loose focusing geometry as well as by much weaker laser pulses, with tight focusing 
geometry, as recently discussed in Refs. [34|[3^ . 

We verified the scalability by simulating HHG in Ar, using a simulation code that includes both laser and 
XUV field propagation effects, see Methods for details. The dipole response was calculated using the Strong Field 
Approximation [^. Figure [^(e) and (f) illustrate HHG using 10fs laser pulses centered at 800nm. Similar as in 
Figure [^(a) and (b), spatiotemporal intensity maps are displayed, showing the evolution of the total field build-up 
along the nonlinear medium for two parameter sets, differing by 77 = 16 (a factor 256 in input energy!). The total 
field above 31.5 eV (i.e. from the 21st harmonic) is represented. It exhibits a train of ultrashort, attosecond pulses. 
The generation parameters led to significant pulse reshaping effects due to plasma formation, implying that the 
generation conditions were not optimized for efficient HHG. The high intensity leads to divergent, ring-like emission 
except at the rising edge of the laser pulse. Again, an almost perfect scaling behavior can be observed, confirming 
Sq —256 

Experimental verification 

To verify the scaling experimentally, we performed pulse compression experiments via filamentation in gases with 
20 fs input pulses (FWHM) centered at 800 nm. The pulse energy was varied in the range of Si = 0.12 — 2.7 mJ and 
spherical mirrors with focal lengths / = 0.5 — 2.5 m were used to focus into an argon-filled tube with a length ap¬ 
proximately twice the respective focal length. We span a pulse energy range of ^25, the highest energy being limited 
by laboratory space constraints. The pulses emerging from the filament were compressed with chirped mirrors and 
fused silica wedges and characterized using the dispersion-scan technique (see Methods for details). Figures [^a) 
and (b) show temporal intensity as well as spectral amplitude and phase for six different input pulse energies. For 
the shortest focal lengths (lowest pulse energy) gas pressure and pulse energy were optimized for maximum spectral 
broadening and good compressibility, while avoiding multiple filamentation. For all other measurement points, focal 
length and gas pressure were adjusted according to the scaling relations, while the pulse energy was used as a free 
parameter to optimize the output spectrum, resulting in input pulse energies very close to the scaling prediction. 
All employed experimental parameters together with fits visualizing the expected scaling trend are displayed in 
Figure [^c). The post-compressed pulse duration as well as the overall characteristics are very similar for all six 
cases, indicating very good scalability of all relevant linear and nonlinear propagation processes within the employed 
parameter range. 

Up (and even down-) scaling HHG has been investigated previously, albeit in a phenomenological way, and in 
many cases without changing consistently all relevant parameters included in our scaling formalism. To make use of 
high input energies, loose focusing geometries have been implemented since the early days. Although it was often 
realized that in these conditions the use of long media (and low pressures) led to higher XUV pulse energies, to 
our knowledge, no rigorous understanding for this experimental observation has been put forward. In the other 
direction, tight focusing geometries, necessary for HHG with laser systems with low input energy (a few tens of 
/iJ), have often been implemented and found to be detrimental for phase matching, and thus for the conversion 
efficiency. A recent experiment using a short medium at high pressure [^, however, shows a similar conversion 
efficiency as with loose focusing, in perfect agreement with our scaling predictions. 
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Figure 3: Experimental filament scaling: a,b, Measured temporal intensity profiles as well as spectral amplitude 
(b, solid lines) and phase (b, dashed lines) for six different parameter sets, shown in c. For better visualization, the 
plotted datasets are vertically off-set from each other. The measurement was performed selecting the broadband 
radiation on the optical axis more than a focal length distance behind the filament. For reference, the input spectrum 
(gray shaded area) is shown in b. The solid lines in c represent fits to the experimental data points, as defined 
by the presented scaling relations, indicating the expected scaling performance for input laser pulse energies within 
and beyond the measured parameter range. The gray data points in c visualize the extrapolated parameters shown 
in Table m 

Discussion and generality of the scaling principle 

We illustrate the scaling possibilities and experimental challenges for the two phenomena discussed in Table 
Starting from typical experimental parameters corresponding to 5in ^ 1 mJ, we apply our scaling relations both for 
filamentation and HHG up to = 500 mJ and, for HHG, down to Cin = 10/iJ and calculate the expected values 
for output pulse energy, gas pressure and focal length. In the case of filamentation, we use the parameters of the 
experiment presented above and assume that 5out = O.l^in for the central compressed part of the filament. In the 
case of HHG, we start from values close to those reported in with a conversion efficiency in Ar equal to 10“^ 
(see also [^). These examples illustrate the feasibility of both up and down-scaling. Even though long geometries 
need to be implemented, few-cycle laser pulses and attosecond XUV pulses with unprecedented energies are within 
reach. Conversely, high gas densities and very tight focusing geometries are required for the efficient generation of 
attosecond pulses at MHz repetition rates (^ . 

Our scaling model does not indicate any limitation for up-scaling. However, for down-scaling, several effects 
leading to deviations from perfect scalability can be identified. First, non-paraxial propagation effects arise at very 
tight focusing geometries (typically at numerical apertures > 0.3). Second, the not perfectly linear dependence 
of K{uj^p) and possibly ^nl on the gas density [see equation (§1 as well as the weak dependence of l/k{uj^p) on 
the density contribute to increasing deviations from perfect scaling (see SI for details). Finally, at high ionization 
levels and high densities, avalanche ionization, a process that critically depends on plasma dynamics and that is 
not scalable according to our model, can set strict limitations (see SI). In extreme conditions, the generated plasma 
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(mJ) 
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1.5 
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15 

I 

Up-scaled 

500 

5000 

0.045 

18.3 

Down-scaled 

O.OI 

O.I 

2300 

0.08 


Table 2: Extrapolation of typical parameters for filamentation and HHG: The experimentally well explored 
regime around Sin ~ I mJ (denoted as typical) is up- and (for HHG) down-scaled. An input pulse length of 20 fs is 
taken into account for filamentation and 40 fs for HHG, with Ar as nonlinear medium. The beam diameter (FWHM) 
before focusing is in both cases ^ 7 mm. 
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can become opaque (for p ^ 70 bar at 800 nm and room temperature, assuming a totally singly-ionized medium) 
and, for processes like HHG, the single atomic response can be affected by the presence of neighboring atoms [4^ . 
However, we estimate that these effects do not play a major role within the parameter ranges typically employed, 
for example, for HHG in gases and for filamentation [see also Fig. [^d)]. 

The presented scaling framework is very general and applies to other processes, involving linear or nonlinear 
electromagnetic wave propagation in gases. The key condition determining if a nonlinear process is scale-invariant 
is the proportionality ^nl oc p. Nonlinear processes which critically depend on plasma dynamics such as avalanche 
ionization or the acceleration of electrons in relativistic light fields are thus not fully scalable according to our 
derived formalism. Furthermore, for processes that make use of the plasma as source of secondary emission, the 
frequency dependence of the secondary radiation upon gas density induces a non-negligible departure from Pnl oc p 
and thus from scale-invariance. Nonlinear interactions that are scalable to a very good approximation include self- 
focusing, self-phase modulation, wave mixing, as well as field ionization, plasma defocusing and processes involving 
stimulated Raman scattering. Similar scaling principles can also be applied for pulse propagation in waveguides 
such as hollow capillaries [^. We expect our results to be of great interest for ultrafast science since we show how to 
extend different nonlinear methods to the new parameter regimes provided by today’s state-of-the-art femtosecond 
laser technology. Our findings are currently being applied to the design of an up-scaled, next-generation attosecond 
source, for the European facility ELI-ALPS (Extreme Light Infrastructure - Attosecond Light Pulse Source). 
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Methods 

Numerical Methods 

Filamentation: For the filamentation simulations, a uni-directional nonlinear envelope equation was numerically 
integrated [^. This type of equation is valid for pulses down to the single-cycle regime [^. Note that the use of 
an envelope equation is not in contradiction to our discussion; the scaling can be formally shown both for envelope 
and field equations. The propagation equation follows from the nonlinear Helmholtz equation, adapting the method 
suggested by Feit and Fleck to introduce uni-directionality [^. Argon is used as nonlinear medium and dispersion 
is treated by a Sellmeier-type expression for the linear refractive index [^. Self-focusing originating from a Kerr- 
type nonlinearity is described by a nonlinear refractive index n 2 = 0.98 x 10“^^ cm^W“^ for Ar at atmospheric 
pressure [^. The interaction with the generated plasma, i.e. absorption and plasma defocussing, is treated in 
terms of a current density, whereas free electrons are generated via field ionization, numerically described within 
the generalized Keldysh-PPT theory (Perelomov-Popov-Terent’ev) (^[4^ . Avalanche ionization is included via a 
Drude-type model with a phenomenological collision time of tc = 190 fs at atmospheric pressure [^. 

HHG: The numerical HHG simulations are performed in two steps: in the first step the generating laser field is 
propagated through the cylindrically symmetric gas medium, taking into account the time- and space-dependent 
neutral and plasma dispersion and optical Kerr effect, as described in Ref. [^. Then, in the second step the 
harmonic field is propagated using the generated high-order polarization as source term calculated using the strong- 
held approximation [^. For the propagation of the harmonic held neutral dispersion and absorption are taken 
into account. The wave equations for both fundamental and harmonic helds are derived in paraxial approximation. 
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In both equations a nonlinear source term proportional to the gas density is assumed. For the fundamental field 
absorption resulting from optical field ionization is taken into account (also scaling linearly with gas density), while 
neutral and plasma absorption are neglected. The numerical code has been verified by comparison with experimental 
results obtained with the intense harmonic beam line at Lund University. 

Experimental methods 

The experiments were performed with a Titanium:Sapphire chirped pulse amplification laser chain, delivering up 
to 5mJ pulses, with 20 fs duration (FWHM), centered around 800 nm, at a repetition rate of IkHz. The beam size 
was adapted to a diameter (1/e^ of the intensity profile) of approx. 11mm with a reflective telescope, followed by 
a 10 mm aperture. The filaments were generated inside an Ar-filled tube using spherical mirrors of different focal 
lengths (/ = 0.5, 0.75,1,1.5, 2, 2.5 m). The tube was terminated with thin (0.5 mm) fused silica windows in order to 
minimize the impact of dispersion and undesired nonlinear interaction. The Ar pressure was set according to the 
scaling law {p = 8,3.5, 2, 0.98, 0.50, 0.32 bar), while the pulse energy was fine-tuned with a A/2-plate followed by a 
transmissive polarizer in order to match the super-continuum spectra originating from the filament to best possible 
agreement for the different combinations of pulse energy, Ar pressure, and focal length. For the lowest input pulse 
energy, a second polarizer plate was installed in order to improve the polarization contrast. Behind the filament 
tube, the bright center of the filament was selected with an aperture and send to a pulse compressor, comprised of 
movable fused silica wedges and chirped mirrors in double-angle configuration (Ultrafast Innovations). The pulses 
were characterized using the dispersion scan technique (d-scan) [^. The d-scan is based on recording the second 
harmonic spectrum as function of dispersion, in our case realized by moving one of the wedges. The spectral phase 
is retrieved with a numerical algorithm that iteratively reproduces the measured d-scan trace (second harmonic 
spectrum vs. wedge insertion) based on a measured fundamental spectrum, see Figure 1 in the SI. For the temporal 
pulse retrievals shown in Figure |^a), a wedge insertion was chosen corresponding to a minimized pulse duration 
and best agreement with pulse profiles retrieved for other parameter combinations. Note that an offset third-order 
phase remains uncompensated in the employed compression scheme, as visible in the d-scan traces, showing a clear 
variation of spectral frequency distribution as a function of wedge insertion. 

Scaling error calculation 

The normalized relative scaling errors for intensity In = In{r,t) and fluence Fn = Fn{r,z) = /^(r, t, z)dt 

displayed in Figure]^ were calculated for the input parameter set {efn^ Pn^ fn) and referenced to the corresponding 
parameter set with four times higher input pulse energy, all other parameters being adjusted according to the scaling 
relations, i.e. {e^n^ ~ Pn-i = Pn/^^ fn-i = 2/n). For the intensity and the fluence the normalized errors can 
be written as: 

drr|J„-/„i| 

. 2/ood^rdrr|F„-F„_i| 

/^d./ydrrF„ + /^d./~drrF„_i- 

The integrals where computed numerically on a finite grid in r and 2 : using a normalized grid size such that 
Ar/Wo = constant and Azjz-^ = constant with Ar and Az denoting the grid unit size. 
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